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Despite of the success of the slave-boson theory in capturing qualitative physics of high- 
temperature superconductors like cuprates, it fails to reproduce the correct temperature-dependent 
behavior of superfiuid density, let alone the independence of the linear temperature term on doping 
in the underdoped regimes of hole-doped cuprate, a common experimental observation in different 
cuprates. It remains puzzling up to now in spite of intensive theoretical efforts. For electron-doped 
case, even qualitative treatment is not reported at present time. Here we revisit these problems 
and provide an alternative superfluid density formulation by using the London relation instead of 
employing the paramagnetic current-current correlation function. The obtained formula, on the one 
hand, provides the correct temperature-dependent behavior of the superfluid density in the whole 
temperature regime, on the other hand, makes the doping dependence of the linear temperature 
term substantially weaken and a possible interpretation for its independence on doping is proposed. 
As an application, electron-doped cuprate is studied, whose result qualitatively agrees with existing 
experiments and successfully explains the origin of d- to anisotropic s-wave transition across the 
optimal doping. Our result remedies some failures of the slave-boson theory as employed to calcu¬ 
late superfluid density in cuprates and may be useful in the understanding of the related physics 
in other strongly correlated systems, e.g. Na 2 ,Co 02 -yH 20 and certain iron-based superconductors 
with dominating local magnetic exchange interaction. 


I. INTRODUCTION 

The slave-particle theory, which splits the physi¬ 
cal electron into various auxiliary bosonic, fermionic 
and even anyonic elementary degree of freedom, has 
been successfully utilized in numerous strongly corre¬ 
lated electron systems, ranged from frustrated quantum 
magnetismfii^ abelian/non-abelian quantum Hall liq¬ 
uid to heavy fermion compounds and high-temperature 
superconductorsAmong them, the most important 
achievement of the U( 1 )/SU( 2 ) slave-boson theory (or 
the ’plain vanilla’ version, the renormalized mean-field 
theory) is that it has at least qualitatively captured 
basic features of the global phase diagram of cuprate 
superconductore.g., the underdoped pseudogap be¬ 
havior, optimally doped strange metal state and over¬ 
doped Landau Fermi liquid state 4 “— 

However, in spite of those extreme success, one of the 
essential physical measurements of any superconductors, 
i.e. the London penetration depth in Meissner effect, has 
not yet been correctly reproduced and explained in terms 
of the mentioned powerful slave-boson techniquei^ Phys¬ 
ically, the temperature-dependent London penetration 
depth/superfluid density directly detects the supercon¬ 
ducting quasiparticle excitation above the highly entan¬ 
gled many-body superconducting ground-state and thus 
can provide unambiguous clue for the gap structure of 
superconducting pairing symmetryJ^ 

The main challenges for the slave-boson theory in ex¬ 
plaining the superfluid density experiments of cuprates 
are that it cannot reproduce the following two ex¬ 
perimental observations, namely, i) the temperature- 
dependent behavior of superfluid density and ii) doping- 


independence of the low temperature linear-T term of 
superfluid density in hole-doped compounds despite of 
the so-called Uemura scaling (superfluid density is pro¬ 
portional to doping level) has been explained by the 
theoryiiirJl Due to these failures, it was believed that 
non-perturbative effects beyond mean-field and Gaussian 
fluctuations should play an essential role even in the well- 
formulated superconducting states. More seriously, since 
almost all calculations in the slave-particle formalism are 
performed in terms of the framework of mean-field theory 
and Gaussian gauge fluctuation, these failures may imply 
the painful breakdown of the slave-particle technique in 
any unconventional superconductivity. 

Furthermore, we notice that in electron-doped cuprate, 
e.g. Nd2-xCea,Gu04 and Pr2_a;Gea;Cu04, it has been 
firmly established that both electron and hole Fermi 
pockets are responsible for the resulting d-wave super¬ 
conducting stateji^^— However, in the framework of slave 
boson theory, even qualitative treatment on their Lon¬ 
don penetration depth is not reported although two 
Fermi pockets behavior. Big, B2g Raman scattering spec¬ 
tra and inelastic neutron scattering spectra have been 
well-explained by slave-boson mean-field theory with as¬ 
sumption of (tt, tt) anti-ferromagnetic spin-density-wave 
(SDW) orderj^^— 

But, considering the fact that the slave-boson theory 
has provided so much interesting physics of cuprates, we 
hold the point that the conventional (Higgs confined) su¬ 
perconducting state should also be well described by this 
theory with appropriate extension and/or modification. 
In this paper, we focus on the mentioned issue of su¬ 
perfiuid density and revisit the formulation of the Lon¬ 
don penetration depth. Instead of the current-current 
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correlation function, we employ the London relation to 
calculate the superfluid density. It is found that the ob¬ 
tained formula can reproduce correctly the temperature- 
dependent behavior of the superfluid density in terms of 
the usual slave-boson mean-field formalism. Such new 
formalism emphasizes the weaken quantum correction 
in the current-current correlation of slave-boson theory 
accompanying with the disappearance of doping depen¬ 
dence in paramagnetic current response part, thus the 
result of the low temperature London penetration depth 
is now consistent with experimental data. More inter¬ 
estingly, at low temperature, the doping dependence of 
the linear temperature term of superfluid density is found 
to be substantially weakened and if a physical supercon¬ 
ducting gap is considered, independence of doping depen¬ 
dence can be readily realized. [Note, however a fluctuat¬ 
ing d-wave superconductor model with weakly interacting 
Bose gas may give rise to the doping- independence;^] 

Moreover, the impurity and multi-band correction are 
analyzed. It shows weak correlation effect on the dy¬ 
namics of impurity scattering. The superfluid den¬ 
sity formula is derived for the electron-doped cuprate 
and the calculation is qualitative consistent with exper¬ 
iment. Furthermore, the approximated low tempera¬ 
ture formula explains the origin of d- to anisotropic s- 
wave transition across the optimal doping, which is a 
success of the present theory. We hope that the find¬ 
ing obtained in the present work may be also relevant 
for unconventional superconductivity in triangular lattice 
compound Naa;Co02-yH20 and certain iron-based su¬ 
perconductors with dominating local magnetic exchange 
interaction;^!^ 

The remainder of this paper is organized as follows. In 
Secini the mean-field Hamiltonian of usual t — J model 
on square lattice is introduced and briefly discussed. In 
secnni the superfluid density formula in the slave-boson 
mean-field theory is reviewed. Next, in Sec lIVI we provide 
an alternative formalism for superfluid density, which 
gives correct result when comparing with experiments. 
In SecjV] we give some discussions on the impurity and 
multi-band effects, which may be important for real ma¬ 
terials like electron-doped cuprate and iron-based super¬ 
conductors. Finally, we end this work with a brief con¬ 
clusion in SecIVIl 

II. THE MEAN-FIELD MODEL 

The standard mean-field Hamiltonian for the t — J-like 
model in the slave-boson framework reads^i^ 

k <7 k 

+2J{x^ +—), ( 1 ) 

where the single-particle energy spectrum including up 
to the third nearest-neighbor hopping is Ek = {—2tS — 
Jx)(cos kx -I- cos ky) — 4t'S cos kx cos ky — 4 t"d(cos^ kx + 


cos^ ky — 1 ) — fi with 6 denoting the doping level. The 
strong correlation effect from prohibiting double occu¬ 
pation on one site in the original t-J model is encoded 
with explicit doping dependence in the energy spectrum. 
The mean-field parameters are defined as {fl^fja) = X 
and (/jt/jq — /zq/jq) “ which are from the decou¬ 

pling of Heisenberg interaction in the particle-hole and 
particle-particle channel, respectively. 

It is important to recall that the definition of these 
mean-field parameters is motivated by the local spin- 
singlet idea of resonance-valence-bond (RVB) quantum 
liquid^, thus the local pairing is formed from the begin¬ 
ning and no super-glue like phonon is involved. How¬ 
ever, the more conventional spin-fluctuation theory can 
also give rise similar effective pairing Hamiltonian EqU 
by extracting singlet paring interaction from the static 
transverse susceptibility, which results from the local 
Heisenberg interaction.'^ Therefore, in this sense, the use 
of slave-boson mean-field formalism does not necessar¬ 
ily link to RVB but only indicates the strong coupling 
feature of the problem itself imposed by the non-double 
occupation condition. 

The physical electron is transformed into Ca = '/Sfn 
with /ct being the fermionic spinon excitation. It should 
be emphasized that although we have used the slave- 
boson representation to rewrite the electron operator 
with fractionalized fermionic spinon, the theory itself 
does not work at deconfined state. In other words, no 
fractionalized quantum phase is involved in our calcula¬ 
tion and all considered states here are ’normal’ phases, 
which could be understood from general Fermi liquid 
framework as enthusiastically advocated by Laughlin;^ 

Then, using the Bogoliubov transformation Ak^ = 
Aifc/fet + = -Vkfkt + Mfe/lfej. and nl = 

1 — v'l = 5 ^I + we can obtain the diagonalized 

Hamiltonian as 

H ^k{Al^Ak-i + A^_i^^A-ki) + “ -®fc) 

k k 

+2Jix^ + ^) ( 2 ) 


with Ek = + Usually, on the square 

lattice, one considers the dx^-y^-wave pairing = 
A(cos kx — cos ky) and the extended s-wave pairing Ak = 
A(cosfca; -I- cosky). For the cuprate, it has been firmly 
established that the dominated pairing symmetry is the 
anisotropic dx^-y^-wave^^ although topological p+ip, 
d+id and even Fulde-Ferrell-Larkin-Ovchinnikov type 
pairing structure might be relevant to the heavily un¬ 
derdoped cuprate;^^— Therefore, we will only consider 
the case of dx2-y2-wave in the remaining part of the 
present paper, and extension to other pairing symmetry 
is straightforward. 
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III. SUPERFLUID DENSITY IN SLAVE-BOSON 
MEAN-FIELD THEORY 


The standard superfluid density formula in the frame¬ 
work of slave-boson mean-field theory reads as follows 
(See Appendix A for details);^ 


m 



+2S 


(dsk^ 

\^dfF{Eu) 

\dkxy 

1 dEk 


( 3 ) 


Here, the first term, which mainly contributes to the zero- 
temperature superfluid density, has negligible depen¬ 
dence of temperature. However, in the second term, the 
thermal excited superconducting quasiparticles deplete 
the superfluid density, which leads to the temperature- 
dependent behavior. It is emphasized that since the sec¬ 
ond term comes from the paramagnetic current-current 
correlation, which involves two single-electron Green 
functions ^ (cc'^)(cc'f) ^ ^^(//^)(//^), the <5^ de¬ 

pendence is expected. However, as what can be seen 
in Fig[I] the normalized superfluid density ns{T)/ns{0) 
versus T/Tc dose not vanish at critical temperature T^, 
which implies the breakdown of the standard superfluid 
density formula. 



FIG. 1. The failure of the superfluid density formulation Eq. 
© obtained by the conventional current-current correlation 
function. The normalized superfiuid density UsiT)/ns{Q) ver¬ 
sus T/Tc at doping level 5 = 0.15 dose not vanish at critical 
temperature Tc. 

More seriously, if we try to inspect the low temperature 
behavior of the superfluid density, it reads^i^ 

(4) 

with ns(0) ^ 5. It has been emphasized by Lee, Na- 
gaosa and Wen that ns(0) ^ 5 agrees with experiments, 
however the doping dependence in the linear-T term ^ 5"^ 
dose not.— As observed in existing experimental data, the 
linear-T term should not show apparent dependence on 
doping level 5^— Importantly, such inconsistency leads 


to the belief that non-perturbative effects beyond mean- 
field and Gaussian fluctuation play an essential role even 
in the well-formulated superconducting states.— Since al¬ 
most all of calculations in the slave-particle framework 
are established in terms of the mean-field theory and 
Gaussian gauge fluctuation, such failure may break down 
even the application of slave-particle technique in any un¬ 
conventional superconductivity. 

Theoretically, one expects fluctuation effect beyond 
mean-field and Gaussian level may cancel out the un¬ 
desirable doping dependence. In this respect, the nu¬ 
merical variational Monte Carlo method, which exactly 
performs Gutzwiller projection on each site for mean- 
field wave-function may give more reliable estimation 
on corresponding physical quantities but calculation on 
temperature-dependent superfluid density has still not 
been reported in literaturei^ 


IV. ALTERNATIVE FORMALISM FOR 
SUPERFLUID DENSITY 


In this section, we will present the details on an alter¬ 
native formalism for superfluid density, which can over¬ 
come the mentioned difficulty. The derivation procedure 
is conventional in fact and can be found in the standard 
text books i^^'^^ The main point is that the no current- 
current correlation is used for this derivation, thus the 
doping effect is weakened. We should remind the reader 
that similar formalism has been successfully used in the 
superfluid density calculation about hole/electron-doped 
cupratei^i^^ and t — J model on the honeycomb lattice— ^ 

First, we define the so-called London relation (The 
-» 2 -* 
scalar case is just J = 

Ji — Qij Aj , 


where Qij is a tensor, which is related to the superfluid 

i 2 

density as Qu = And, Ji and Ai represent the 

electronic current density and external electromagnetic 
vector potential, respectively. Meanwhile, the current Ji 
has two components, namely, the paramagnetic part J* 
and the diamagnetic part J^. 


A. Paramagnetic current J* 


By the slave-boson technique, the paramagnetic part 
j; = -nev = -eEfcK(4tCfet)+^^-fe(c-fe;C_fc^)) can be 
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written as 


k 

= -e6j2vli{AlA,^) - {Al,^A.,^)) 

k 

=-eSY^vKlFiEk) - fF{E-k)) 


= 2e^^E 

k 



dfF{Ek) \ 
dEk j 


Ar 


In the above derivation, we have used the relation that 
when there is the external vector potential A, fFiEk) — 
friE-k) ^ -2e^^^^vlAj. Thus, we have = 
2e^5 X]fc ['^fc^fc ] • (The effective velocity is defined 

by vl. = §^-) It is clear to see that in contrast to the 
conventional formula in last section, the dependence on 
doping level is ^ (5 rather than 5^, thus the encoun¬ 
tered doping dependence problem actually disappears. 
The reason of this distinction is that here the contribu¬ 
tion of the paramagnetic current to the superfluid re¬ 
sponse is effectively single-particle-like rather than the 
two-particle correlation in the usual current-current cor¬ 
relation function. Usually, as seen from the more sophis¬ 
ticated variational Monte Carlo calculation^ the slave- 
boson mean-field theory often overestimates the effect of 
strongly electronic correlations and the result presented 
here suppresses partial unphysical degree of freedom in¬ 
troduced in slave-particle representation. 


B. Diamagnetic current 


Correspondingly, the diamagnetic part = 


+ (clkiC-ki)) 

‘■fc 


1 

can be written as 


J^= ^ ^ [{fl^fk,) + ifUf-H)) 


= -e^6A,J2 


— {f^k - ^k){{A\^Ak^) A {A^_f.^A-k])) 


E^2ut 


= -e^SAj 

k 

= -e^5Aj 


mi 


1 £fc 

m''^ Ek 


VF{Ek) 


mi 


^k 


_ 1_ EEtanh(^) 

11 \ T7' ^ / 


E(i - 

71*^ Ek 


Ek 


E, 


2T' 


Therefore, the diamagnetic kernel reads Qfj = 
(l- ^tanh(^))] 


C. The superfluid density n® 


Since the total superfluid response kernel Qij = + 

Qy , we have 



£^tanh(:^)) 
Ek 2T’’ 




dfFjEk) 

dEk 


and 


Q^i = e^6 


k 




dfFjEk) 

dEk 


So, the superfluid density reads 

= ^E 


rri 




d'^e 


Ajvl) 


2dfF{Ek) 


dEk 


fe/1 Efe 


dkf Ek 2T^^ ^ dh ^ dEk 


Usually, since the considered system is symmetric be¬ 
tween X and y direction, which appears in the slave-boson 
mean-field theory on square lattice, we can use the sim¬ 
plified formula below 


m 



l_£itanh(^) 
Ek ^2T> 


n(dek.2 dfF{Ek) 
^dh’ dEk 


(5) 

Obviously, one can observe that the undesirable doping 
dependence does not appear in the quasiparticles deple¬ 
tion term (the second term), which leads to the phys¬ 
ically correct result as shown in Fig. and qualita¬ 
tively agrees with the penetration depth measurement in 
hole-doped cuprate superconductor Furthermore, 

the formulation obtained can be employed to all dop¬ 
ing cases, ranged from underdoped to overdoped regime, 
as shown in Fig. [3] [Obviously, since no peusdo-gap is 
introduced in the present model, the comparison to un¬ 
derdoped cuprate should not be considered seriously. It 
seems that the phenomenological theory in Refl^ may 
be useful in treating the issue of peusdo-gap.] The result 
is also consistent with reported data in Refs. [3 and 
[^ . In other words, the curvature of the normalized 

superfluid density ^ ( nl[oj ) nearly identical and is 
independent of doping level, thus recovers the observed 
behavior in experiments. 

Consider the low temperature behavior of the super¬ 
fluid density, we have 


ns(T) Ri ns(0)(l 


(2ln2)r 

A 


)■ 


( 6 ) 
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FIG. 2. The qualitatively correct temperature-dependence 
of the superfluid density given by Eq. The normalized 

superfiuid density ns{T)/ns{0) versus T/T^ at doping level 
S = 0.15 by using obtained superfluid density formalism. 






FIG. 3. The normalized superfluid density na{T)/ns{0) at 
different doping levels ranged from underdoped to overdoped 
regimes. 

Comparing with Eq. 0, the doping dependence of linear 
term is weakened to ^ 5 instead of ^ <5^, which improves 
significantly the slave-boson theory at low temperature. 
However, both the experimental data and phenomeno¬ 
logical theoretical analysis of superfluid density suggest 
the following low temperature formula)^^— 

n,(r)-n.(0)^-^, (7) 

which shows no apparent dependence of doping level. 
(The gap A(0) may only have weak dependence of dop¬ 
ing.) If one considers the fact that in the slave-boson 
theory, the true superconducting gap of original electron 
is Asc = SA, one has 

ns{T) - ns{0) « ^ - J. (8) 

This formula is identical to Eq. ©. Therefore, an im¬ 
proved slave-boson theory description with the physical 


pairing gap is able to reconcile the difficulty between the 
present theory and the experimental measurements. 

V. IMPURITY AND MULTI-BAND EFFECT 

For realistic materials, the nonmagnetic impurity scat¬ 
tering is ubiquitous and plays a crucial role in the low 
temperature behaviors of any superconductors, partic¬ 
ularly for nodal pairing states like the most important 
d-wave state. On the other hand, in many real-life mate¬ 
rials, many energy bands contribute to the ultimate su¬ 
perconductivity, thus it is also helpful to investigate the 
formalism of superfluid density in the multi-band case. 

A. Non-magnetic impurity 

Following Ref. , when considering the effect of non¬ 
magnetic impurity on the nodal d-wave state, the corre¬ 
sponding low temperature superfluid density reads 

(9) 

where the zero-temperature superfluid density ns(0) ~ 
n,(0)/(l + i.ssyrVA) is suppressed by the impurity 
scattering T and there is a crossover temperature T* Ri 
O.SSVTA. Furthermore, if T <C T*, the behavior can 
be seen in the superfluid density instead of the linear be¬ 
havior at elevated T T*. It is noted that the only 
effect of strong electron correlation is ns(0) ~ 5. This is 
because that the dynamics of impurity scattering is fully 
encoded by the auxiliary fermions without involving any 
real charge fluctuation. 


B. Multi-band effect: Application to 
electron-doped cnprate 

Generally, the hole-doped cuprate permits a single¬ 
band description and the formalism presented in the main 
text is applicable. 

For the case of electron-doped cuprate, e.g. 
Nd 2 _ 2 :Cea;Cu 04 and Pr 2 _a;Cea;Cu 04 , it has been 
firmly established that both electron and hole Fermi 
pockets are responsible for the resulting d-wave super¬ 
conducting statei^ Theoretically, the t — J model has 
been utilized with inverse sign of hopping parameters 
as compared to their hole-doped counterpart!^ The 
observed two Fermi pockets behavior is captured by 
the slave-boson mean-field theory with with assumption 
of (tt, tt) anti-ferromagnetic spin-density-wave (SDW) 
order Other physical quantities like Big, B 2 g Raman 
scattering spectra and inelastic neutron scattering 
spectra have also been well-explained by slave boson 
theory i^^i^'^ However, we emphasize that there is still 
no qualitative treatment on the London penetration 
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FIG. 4. The normalized superfluid density ns(T)/usiQ) at 
optimal doping levels 5 = 0.13 for electron-doped case. 


depth in those electron-doped cuprate. [Note, however a 
phenomenological two-band model has succeed in fitting 
and explaining on the experimental data by our previous 
workJiSi] 

Here, utilize the t — J-like model in previous works and 
extend the discussion in last section,i22r21 we can obtain 
the following superfluid density formula. 





dki 




k \2 


dfriE^) 

dE- 


Here, = {sk + Efc+g ± \/{sk - efe-rg)^ + {Jm)^)/2 de¬ 
notes two antiferromagnetic energy bands, which results 
from the preformed anti-ferromagnetic SDW long-ranged 
order with Q = (7r,7r)j^ If no doping is introduced , the 
system is an anti-ferromagnetic Mott insulator. Upon 
doping, the band gives rise to the electron-like Fermi 
surface centered at (7r,0). When approaching optimal 
doping level, the band drives the formation of hole-like 
Fermi surface around (7r/2,7r/2) and the superconducting 
instability ultimately develops with appearance of such 
hole Fermi surface. The correspon ding superco nducting 
quasi-particle spectrum is E^ = yj with the 
assumption of d^ 2 _y 2 -wa,ve gap Ak^ In FigjU we have 
shown the normalized superfluid density at doping levels 
6 = 0.13 (optimal doped), which is consistent with the 
published data of the penetration depth measurement in 
Refill. 

More importantly, when inspecting the low tempera¬ 
ture behavior of the superfluid density, we find 

MT) = <(0)(1 - + n;(0)(l - h^). 

Here, A' is the minimum value of gap function Ak at 
the (tt, 0) electron Fermi surface with constant a and b. 
With this approximated formula, one can see that when 
doping is not large (n+(0) nj(0)), only the first term 

dominates, which leads to the gapped s-like exponen¬ 
tial behavior observed in experiments When the hole 
Fermi surface forms around optimal doping, the latter 
term with linear T behavior competes with the first term, 
thus behaves like a gapless nodal superconductor. There¬ 
fore, this simplified formalism, which agrees well with the 
results in our previous studysuggests that there is no 
true d- to anisotropic s-wave transition across the opti¬ 
mal doping 


C. Multi-band effect: Iron-based superconductors 

When entering the ‘iron age’, it has been suggested 
that certain iron-based superconductors can be viewed 
from the perspective of the doped Mott or orbital selec¬ 
tive Mott insulator, whose pairing originates from the lo¬ 
cal magnetic exchange interaction in the parent FeAs and 
FeTe compoundsThus, their pairing physics may 
be explored by multi-band version of t — J model since 
more than one Fe d-orbital contribute to the low energy 
physics. Particularly, the superfluid density formula de¬ 
rived in the main text can be used in this case with proper 
extension to include orbital effect. 


VI. CONCLUSION AND DISCUSSION 


We have provided an alternative formalism for super¬ 
fluid density, which can reproduce qualitatively the exist¬ 
ing experimental data in high-T^ hole and electron-doped 
cuprate superconductors. This new formalism resolves 
the failure of slave-boson theory as applied to calculate 
the temperature-dependent superfluid density in t — J- 
like strong coupling lattice fermion models. It also indi¬ 
cates that the slave-boson mean-field theory generically 
overestimates the many-particle correlation effect in the 
many-point correlation function, which should be further 
studied by more sophisticated techniques like variational 
Monte Carlo simulation and cluster dynamic mean-field 
theory. Therefore, although only superfluid density is an¬ 
alyzed in the present work, we believe that other physical 
observable dominated by many-particle correlation, e.g. 
optical conductance and Hall coefficient, should be cal¬ 
culated with similar formalism as ours when slave-boson 
theory is applied. 
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Appendix A: The standard snperflnid density 
formula in slave-boson mean-fleld theory 

The standard superfluid density n® formula reads as 
follows!^ 

5 ^ ^ (A1) 
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where (•) means the expectation in the BCS mean-field 
Hamiltonian and 


ka 

k 

= Sj}{r). (A2) 


Thus, we have 


= J dT{j}{0)j}{T)). 

(AS) 

After calculating the expectation in the above equation, 
the final result reads 



TO* 






dek^2 


d/pjEk) 

dEk 
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